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We consider black hole production at the LHC in a generic scenario with many extra dimensions 
where the Standard Model fields are confined to a brane. With ~ 20 dimensions the hierarchy 
problem is shown to be natnrally solved without the need for large compactification radii. We find 
that in such a scenario the properties of black holes can be used to determine the number of extra 
dimensions, n. In particular, we demonstrate that measurements of the decay distributions of such 
black holes at the LHC can determine if n is significantly larger than 6 or 7 with high confidence, 
and thus can probe one of the critical properties of string theory compactifications. 


One of the most difficult questions facing theoretical 
high-energy physics is how to consistently combine Gen¬ 
eral Relativity with Quantum Mechanics, as naive quan¬ 
tization produces unrenormalizable divergences. This is¬ 
sue is exacerbated by the hierarchy problem, which asks 
why the electroweak scale, ~ TeV, is so small com¬ 
pared with the (reduced) Planck scale, Mpi ^ a, few 10^® 
GeV, which is associated with the energy at which non- 
renormalizable Einstein gravity becomes strong. It ap¬ 
pears that resolution of these puzzles may require a com¬ 
plete theory of quantum gravity. 

In recent years it has been proposed that the funda¬ 
mental scale of gravity might not be Mpi, but rather 
M* ~ TeV P, U- There is then no large hierarchy be¬ 
tween the gravitational and electroweak scales. In this 
scenario, the observed weakness of gravity results from 
the presence of extra dimensions with large radii. In the 
simplest picture, gravity is able to propagate in all D 
dimensions, but the Standard Model (SM) fields are re¬ 
stricted to a 3 -I- 1 dimensional “brane”. The strength of 
gravity at long distances is then diluted by the volume of 
the extra dimensions. Here, we examine a scenario where 
the number of extra dimensions is large. In this case, as 
we will see below, additional hierarchies do not arise be¬ 
tween M* and the size of the additional dimensions. In 
particular, we examine the properties of black hole (BH) 
production and decay at the LHG with different num¬ 
bers of extra dimensions and show that the number of 
additional dimensions n can be determined at high con¬ 
fidence, in particular when n is large. Our results hold 
in the generic case where the size of the BH is much 
less than the curvature of the additional dimensions, and 
where the SM is confined to a 3-brane. 

As of now, the best candidate for a complete theory of 
quantum gravity is (critical) string theory (GST), which 
reduces to Einstein gravity at low energies and allows 
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for the computation of finite S'—matrix amplitudes. For 
GST to be a consistent theory there are three essential 
ingredients: (i) the fundamental objects of the theory 
are no longer point-like and must have a finite size of or¬ 
der Mg, the string scale; (ii) supersymmetry must be a 
good symmetry, at least at scales > Mg] {Hi) space-time 
must be ten or eleven dimensional, {i.e., D = A + n = 10, 
if the string coupling is perturbative, D = II if it is 
non-perturbative), with the additional dimensions being 
compactified at a radius Rc ^ A/Mg. Most research in 
string theory so far has focused on critical string theo¬ 
ries, where the world-sheet anomalies are automatically 
canceled. It is precisely this anomaly cancelation that 
requires D = 10. However, there are consistent non- 
critical backgrounds of string theory in arbitrary num¬ 
bers of dimensions. Here, the anomalies are canceled by 
solving the equations of motion taking into account the 
tree level moduli potential as well as contributions to the 
equations of motion from other sources such as fluxes, 
orientifolds, and branes |^. In either case, the common 
expectation is that Mg is slightly below or equal to Mpi 
which would imply that the predictions of GST are dif¬ 
ficult to test directly. Currently there is no evidence for 
any of these basic assumptions. If indeed Mg ^ Mpi it 
may be that GST can never be directly tested in labora¬ 
tory experiments. Furthermore, even if supersymmetry 
and/or extra dimensions were discovered in future exper¬ 
iments, this would be no guarantee that CST represents 
the correct theory of nature. 

We will show in this paper that the number of com¬ 
pactified large dimensions can be determined from black 
hole production at the LHC. This would provide a probe 
of classes of CST models. Specifically, if n > 6(7) is 
measured with high confidence then present CST com- 
pactffications would be tested. As a proof of principle 
for our proposal, we will show that there exists a region 
in the parameter space where we can experimentally ex¬ 
clude the case n < 6(7) at 5a significance. 

For purposes of demonstration, we perform our cal¬ 
culations in the the large extra dimensions picture of 
Arkani-Hamed, Dimopoulos and Dvali (ADD) [^. We 
emphasize that our results are general and we only use 
ADD as a calculational framework. Here, M* and Mpi 
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are related by Mpi = VnM^^'^, where 14 is the volume 
of the n compactified large dimensions. For simplicity 
in what follows, we will assume that this n—dimensional 
space is compactified on a torus of equal radii so that 
14 = (27ri?c)”, where Rc is the compactification radius. 
Given Mpi and M* ^ a few TeV, Rc becomes completely 
fixed by the relation above. Note that the case n = 1 is 
excluded while n = 2 with low M* is disfavored by cur¬ 
rent data 0. For the case of a torus, the graviton has 
Kaluza-Klein(KK) excitations , with masses given by 
where n labels a set of occupation num¬ 
bers. The KK graviton couplings to the Standard Model 
(SM) fields are described by the stress-energy tensor 
given in D dimensions by £ = — 

This scenario has three distinct experimental signatures 
which have been studied in some detail in the literature: 
(i) missing energy events associated with KK graviton 
emission in the collisions of SM fields; {ii) new contact in¬ 
teractions associated with spin-2 KK exchanges between 
SM fields {in) black hole production in particle col¬ 
lisions 

Is there any guide as to what values of n > 6(7) we 
should consider? For n < 6 it is well known that the 
hierarchy problem is not truly solved. Although we have 
reduced M* to a few TeV, M^:Rc » 1, as seen in Fig 
n By contrast, with n large we could have M^,Rc ^ 10. 
Note that, if M^Rc < 1 the theory would lose its pre¬ 
dictive power since the compactification scale is above 
the cutoff. To obtain the interesting range of compact¬ 
ification radii, 1 < M^Rc ^ 10, requires 17 < n < 39, 
hence we will focus on this set of values in what follows. 
If the compactification topology is a sphere, rather than 
a torus, this changes to n ^ 30, as seen in Fig^ It is 
important to notice that this model does solve the hi¬ 
erarchy problem for large n, but this would lie outside 
the realm of CST. Note that some other modifications 
of the compactification geometry can obtain < 10 

[m For such large values of n the Kaluza-Klein masses 
are at the TeV scale. For example, in the ADD case, 
since each graviton KK state is coupled with 4 dimen¬ 
sional Planck strength, Mpi, to the SM fields, it is clear 
that this significantly weakens the KK contributions to 
the processes {i) and (ii) above. Thus, no meaningful 
collider constraints would be obtainable; this may also 
happen in the generic model we consider here. For ex¬ 
ample, in ADD with n = 2, precision measurements at 
the International Linear Collider at = 1 TeV will be 
sensitive to M* < 10 TeV, while with n = 21, this drops 
to M* ^ 1 TeV. Thus for reasonable values of M* the 
only signal for large n in ADD is black hole production. 

We now investigate BH production at the LHC in de¬ 
tail; for previous studies see na. When > M* BHs 
are produced with a geometric (subprocess) cross section, 
a ~ 7ri?g. We expect this to hold in all models which sat¬ 
isfy our assumptions. Here Rs is the Schwarzschild radius 
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FIG. 1: MfRc as a function of n for M* = 1 TeV for a torus 
(solid) and sphere (dashed) compactifications. 


corresponding to a BH of mass MpH — VS- Rs is given 

by 0 


M^Rs 
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( 1 ) 


Note that a n for large n. Numerical simulations and 
detailed arguments have shown that the geometric cross 
section estimate is good to within factors of a few [l3| . 
The total number of BH events at the LHC with invari¬ 
ant mass above an arbitrary value Mbh min is shown in 
Fig.121 The scale of the total inclusive BH cross-section, 
~ 100 pb, is huge compared to that which is typical of 
new physics processes, < 1 pb. Thus, over much of the 
parameter space the LHC will be producing over a mil¬ 
lion BH events per year. This high rate means that there 
will be tremendous statistical power, and essentially all 
measurements will be systematics limited. 

The semiclassical treatment, used here and in all previ¬ 
ous studies 0 , may recieve potentially large corrections 
from two sources: (f) distortions from the finite com¬ 
pactification scale as Rg approaches Rc, and {ii) quan¬ 
tum gravity. Case (?) is easily controllable. We know 
that in 5 dimensions the critical point for instabilities 
due to finite compactification is {Rg/Rc)^ ~ 0.1 0 . For 
LHC energies we always have {Rg/Rc)"^ <C 0.1, so these 
corrections are negligible. In more dimensions the ratio 
of the volume of a BH with fixed Rg to the volume of 
the torus with fixed Rc drops rapidly with n, so we ex¬ 
pect the corrections to be even smaller. Case {ii) is more 
problematic; we estimate the quantum gravity effects by 
looking at the corrections from higher curvature terms in 
the action, e.g. 


S = 


M. 


D-2 






• ( 2 ) 


Here R is the Ricci scalar, and Ci is the ?th order Love¬ 
lock invariant, with £2 being the Gauss-Bonnet term |l6| . 
This equation also defines our convention for the funda- 
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FIG. 2: Top panel: Cross-section for production of black 
holes with mass M > with M* = 1.5 TeV, for 

n = 2(bottom) to 25 (top) of the band. Also shown is the QCD 
dijet cross-section for dijet invariant mass M > and 

I 77 I < 1. Bottom panel: distribution of BH events passing 

cuts described in the text for M* = 1 TeV and n = 2, 6 , 21. 


mental scale Schwartzchild solutions are known 

for arbitrary values of the ai 13 If we assume that 
the higher curvature terms are radiatively generated, and 
hence each ai is the ith power of an expansion parameter 
a (as occurs in string models 0 ) , we find that aD’^ < 1. 
For a of this size we find that the corrections are always 
less severe as n increases, with a ~ 20% correction to Rs 
for n = 20. This does not qualitatively affect our conclu¬ 
sions here; for a more detailed study of these corrections, 
see 0113 . 

We now come to the crucial question, is there any prop¬ 
erty of the produced black holes that can resolve the num¬ 
ber of dimensions? The cross-section is n-dependent, but 
the overall scale is set by 1/M^, so one would first have 
to measure M* independently to good accuracy to ob¬ 
tain any resolution on n. Cross section ratios at different 
BH masses could be used, however, the range of energies 
that are clearly in the geometric regime and accessible 
to the LHC is not likely to be large. This leads us to 
the decay properties of black holes. One generically ex¬ 
pects that black holes produced at colliders are formed in 
highly asymmetric states, with high angular momentum, 


^ We note that this is related to the other definitions in the lit¬ 
erature by M* = (Stt) '•+2Mdl|3 = [ 2 ( 2 ' 7 r )0 ^’•+ 2 A 1 'gtE 3 
= (27r)“"/(2-l-")yf^^ as used in Giudice et. al |^. 


and possibly a non-zero charge. However, they quickly 
shed their charge and angular momentum by emitting 
bulk graviton modes and soft brane modes, and relax to 
a simple Schwartzchild state; their decay then proceeds 
primarily by thermal emission of Hawking radiation 0 
until Mbh where quantum gravity effects will me¬ 

diate the final decay. The Hawking temperature is given 
by 


Th = 


{ji -\- 

47r 


r(n^) 


Mbh 


{n -h 2)7r(”+3)/2 M* 


-l/(n+l) 


( 3 ) 


From this we can see that, at fixed Mbh, higher dimen¬ 
sional BHs are hotter. Since the average multiplicity goes 
inversely with the temperature, a low dimensional BH 
will emit many quanta before losing all of it’s energy. 
By contrast, the decay of a high dimension BH will have 
fewer final state particles, and each emitted quanta will 
carry a larger fraction of the BH energy. We will use 
this difference to obtain experimental resolution on n. It 
was seen in 0 that for n < 6 an error of ±0.75 could 
be obtained. However, as n gets large the BH properties 
at adjacent n converge, so it is a-priori unclear at what 
level n can be determined, if at all, in this case. 

The previous argument suggests we examine the final 
state multiplicity, or the individual particle px distribu¬ 
tions as a probe of n. The multiplicity is affected by two 
major sources of uncertainty: (a) contributions from ini¬ 
tial and final state radiation that produce additional jets, 
and (b) the details of the final quantum gravity decay of 
the BH are unknown. In what follows we will assume that 
this remnant decay is primarily 2-body. However, this is 
clearly model-dependent; we prefer observables that are 
independent of this assumption, disfavoring the multi¬ 
plicity. By contrast, the px spectra of individual parti¬ 
cles, particularly at high-pT, will be mostly sensitive to 
the initial temperature of the BHs. There are many such 
distributions that one could consider. In particular, one 
would like to examine all possible distributions and see 
that the candidate BH states are coupling equally to each 
SM degree of freedom, verifying that these are gravita¬ 
tional phenomena [3 • For illustration we will focus here 
on the ]^x and individual jet px distributions for the BH 
final state. 

To calculate these distributions, we have simulated BH 
events using a modified version of CHARYBDIS [^ . 
linked to PYTHIA First, a large sample of BHs 

with masses above a critical value Mmin = M^ is gen¬ 
erated. From these we select events by cutting on the 
reconstructed invariant mass, Mjnv of the event, defined 
by summing over all visible final state particles or jets 
with rapidity \r]\ < 3, and with px > 50 GeV. We would 
like to select events where the BH mass is large enough 
that the event is in the geometrical regime, and quantum 
gravity corrections are small. To do this, one would need 
to extract from the data an estimate of the size of M*. 
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While we have no fundamental model for the quantum 
gravity effects near threshold, we can assume that there 
will be a turn-on for BH production near M*, and the 
cross-section will then asymptote to the geometric value. 
While this will not lead to a precision determination of 
M*, it can clearly be used to set an optimum cut on 
Minv In the context of a particular model of the thresh¬ 
old based on the action we find that Minv > 2M* is a 
reasonable cut M- We include initial-state radiation in 
the simulations, since that can lead to a contamination 
of lower \/l events in our sample. In the case of jets, for 
simplicity we turn off hadronization, and simply look at 
the parton-level characteristics. 

To be specific, we generate a “data” set of ^ 300A: 
events with n = 21 and M* = 1 TeV. We use this size 
sample as a conservative lower estimate of BH produc¬ 
tion. If the cross section is within an order of magnitude 
of that in Fig. |21 the LHC will collect many millions of 
events, giving an increase in statistical power over that 
presented here. Alternatively, if we employed a stiffer cut 
on the lower value of Minv , this would yield a lower statis¬ 
tical sample, similar to the size of 300fc events considered 
here, and we would expect our results to then qualita¬ 
tively hold in this case as well. These “data” events are 
then compared to a number of template sets of events. 
We then ask at what confidence the template can be ex¬ 
cluded by performing a test using only the resulting 
distribution (shown in Fig. We examine the range 
2 < n < 21, and 0.75 < M* < 5 TeV. The lower bound 
on M* comes from non-observation at the Tevatron and 
cosmic rays |23j , while the upper bound is set by demand¬ 
ing that the LHC be able to collect at least 50fc events 
given the cross-section uncertainties. We then determine 
whether the CST region can be probed at high confidence 
within this scenario. For this test case, we find at least 
a 5(7 exclusion for the entire CST region using the pT 
distribution alone, or ~ 40(T using the jet-pr spectrum. 
Though the statistical power in jets is much higher, it 
suffers from more systematic uncertainties. Fig. 13 shows 
the 3, 5, and lOa exclusion contours in the (n, M^,) plane 
obtained using the distribution for this test case. If 
the LHC collects a few million events rather than the 
300fc sample used here, simple scaling tells us that the 
5a curve excludes n < 20, and the IOct curve excludes 
n < 11. 

We have shown that the CST region can be excluded 
within this scenario if n = 21. What about other values 
of nl On changing the number of dimensions used in 
generating the “data”, we find that for any n > 15 the 
CST region can be excluded by at least 5tT, with 300fc 
events. We would, of course, like to know in what region 
of the parameter space this type of definitive test can 
be performed. A more detailed study of the parameter 
space is in progress [T^. 

In conclusion, we have shown that if there exist 
many TeV sized extra dimensions with the SM fields con- 
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FIG. 3: Exclusion curves in the (n,M*) plane, assuming the 
data lies at the point (21,1 TeV). Points outside the curves 
are excluded at 3, 5, or 10a. 


fined to a 3-brane, then there exists an observable that 
can probe classes of critical string theory models. 
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